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Si(lll) and Ge(lll) surface formation energies were calculated using density functional theory 
for various biaxial strain states ranging from —0.04 to 0.04, and for a wide set of experimentally 
observed surface reconstructions: 3x3,5x5, 7x7 dimer-adatom-stacking fault reconstructions and 
c(2 x 8), 2 x 2 and \/3 x v3 adatoms based surfaces. The calculations are compared with scanning 
tunneling microscopy data obtained on stepped Si(lll) surfaces and on Ge islands grown on a Si(lll) 
substrate. It is shown that the surface structure transformations observed in these strained systems 
are accounted for by a phase diagram that relates the equilibrium surface structure to the applied 
strain. The calculated formation energy of the unstrained Si(lll)-9 x 9 dimer-adatom-stacking fault 
surface is reported for the first time and it is higher than corresponding energies of Si(lll)-5 x 5 
and Si(lll)-7 x 7 dimer-adatom-stacking fault surfaces as expected. We predict that the Si(lll) 
surface should adopt a c(2 x 8) reconstruction when tensile strain is above 0.03. 



I. INTRODUCTION 

During heteroepitaxial growth, elastic strain appears 
when the substrate and growing film lattice constants 
mismatch. This is evident during Ge/Si growth and in 
the vicinity of surface defects like step edges and other 
sites where preferential relaxation occurs. Strain also has 
great impact on surface structure, diffusion of adsorbed 
atoms, thin film growth, nanostructure formation and 
stabilization of facet planes. 1 

A typical example of a strained system is germanium 
on silicon. Interest in Ge heteroepitaxial growth on Si 
stems from several reasons. First, it has shown great 
potential for new high frequency semiconducting de- 
vices, while keeping compatibility with the mature Si- 
based complementary metal-oxide-semiconductor tech- 
nology. Second, the recent success in growing struc- 
turally perfect Ge/Si nanowires and nanorings, may pro- 
vide the structural building blocks for future nanoelec- 
tronic devices. 2 Third, it is a prototypical system for 
studying strained Stranski-Krastanov (SK) growth. Ac- 
cordingly, after the first few smoothly grown epitaxial 
layers (wetting layer), growth proceeds with the forma- 
tion of islands (SK growth mode). 3 This is accompanied 
by several reconstruction changes on top of the growing 
Ge islands. 4 However, technological difficulties exist due 
to the fact that Si and Ge lattice constants differ by about 
4%. This leads to a build-up compressive stress on Ge, 
making the growth of defect-free Ge films on Si a rather 
difficult task. 

Although it is commonly accepted that an applied elas- 
tic strain can induce surface structure transformations, 
precise data quantifying this effect are yet unavailable. 5 
The reported calculated data for Si(lll) and Ge(lll) are 
either based on simplified surface reconstruction mod- 
els (which did not include adatoms and other important 



structural features) or they were obtained using empirical 
potential and tight binding calculations. 6,7 These meth- 
ods not only are limited with respect to their prediction 
accuracy of surface absolute formation energies, but also 
can hardly account for the strain response of complex sur- 
face reconstructions and hybridizations. Hence, the need 
for further progress in this field, in particular the lack of 
surface-strain phase diagrams for Si(lll) and Ge(lll), 
calls for a study of these effects by means of state-of-the- 
art first principles calculations. 

Having said that, this work aimed at establishing a cor- 
respondence between Ge(lll) and Si(lll) surface struc- 
ture changes and the applied lateral strain by combining 
density functional theoretical (DFT) calculations with 
scanning tunneling microscopy (STM) data. Si and Ge 
surfaces were calculated independently. With this ap- 
proach we overlooked any species intermixing that may 
occur at the interface during growth. Nevertheless, given 
the chemical resemblance between both species, this ef- 
fect should have a limited impact. 

The article is organized in the following manner - in 
Section II we describe the computational details, con- 
vergence tests to the calculations, as well as the experi- 
mental procedure to obtain the reported STM data; Sec- 
tion III reports the calculations of surface energies and 
surface strain response, as well as a discussion of these 
results with the STM data as scenery; and finally we 
summarize the conclusions in Section IV. 



II. METHODS 
A. Computational details 

Absolute surface formation energies were calculated ac- 
cording to the recipe by Stekolnokov et al. 8 , here adapted 
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for strained surfaces. Accordingly, the surface formation 
energy of some reference surface should be calculated us- 
ing a slab with two equivalent surfaces (symmetric slab) . 
Relaxed unreconstructed lxl surfaces of Si(lll) and 
Ge(lll) were chosen for such references. The layer in the 
middle of the slab was kept frozen, while atoms in other 
layers were allowed to move without any constraints dur- 
ing atomic optimizations. The surface formation energy 
per unit area (7) is obtained from the energy excess of 
the slab compared to that of a bulk calculation with same 
number of atoms. Hence, for reference surfaces 



Tixi(e) = 
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25ixi(e) 
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where fj, is the Si or Ge strain-dependent chemical po- 
tential, i.e. the energy per atom in bulk under strain e, 
5ixi is the area of a 1 x 1 surface unit cell and -E^ot 
is the total energy of the slab comprising N atoms per 
simulation cell. The factor of 1/2 avoids double counting 
both slab surface energies. 

According to Ref. 8, the calculation of the energy 
gain due to surface reconstructions requires using hydro- 
genated slabs, where one surface layer (hereafter referred 
as bottom layer) is saturated with hydrogen atoms, while 
the reconstructions are considered at the opposite sur- 
face. Two slabs should be used in this approach - one 
for the chosen reference (in our case the unreconstructed 
lxl surface), and a second slab for the reconstructed 
surface under scrutiny. With this setup the location of 
H and Si (Ge) atoms at the bottom layer is kept frozen 
during atomic optimizations, while all other atoms are 
freely allowed to relax. Thus, the energy gain per unit 
area is now 
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FIG. 1. (Color online) (a) Surface formation energy as a func- 
tion of the number of Si(lll) layers in symmetric slabs with 
lxl and 3x3 DAS surface reconstructions, (b) Si(lll)-3 x 3 
DAS surface formation energy as a function of applied biaxial 
strain eb, calculated using a 6-layer thick hydrogenated slabs 
with DZP and mixed SZ-DZP basis sets. 



A 7 rcc( e ) = ~ El^(e)M - n{e)K] , (2) 

where S rec is the unit cell area of the reconstructed slab, 
M = SVcc/Sixi is the number of 1 x 1 reference cells per 
reconstructed unit cell, and K accounts for the number of 
Si (Ge) surface atoms in excess to those in the reference 
cell. The absolute surface formation energy of a recon- 
structed surface under strain is obtained after combining 
the energy reference from the unreconstructed lxl sur- 
face with the energy gain due to reconstructions: 

7rcc(e) = 7ixi(e) + A7 rec (e). (3) 

Total energies were calculated from first principles by 
using the local density functional siesta code. 9 The ex- 
change and correlation functional is that of Ceperley- 
Alder as parametrized by Perdew and Zunger. 10 The 
k-space integrals over three- or two-dimensional Bril- 
louin zones (BZ's) were approximated by sums over 
Monkhorst-Pack (MP) type special points. 11 Norm- 
conserving pseudopotentials were employed to account 



for electronic core states, 12 whereas valence states are 
spanned with help of linear combinations of numerical 
atomic orbitals of the Sankey-Niklewski type, generalized 
to be arbitrarily complete with the inclusion of multiple 
£ orbitals and polarization states. 9 For convergence con- 
trol purposes, large basis calculations were carried out 
by assigning polarized double-^ functions (DZP) to all 
atoms in the slab. This means two sets of s and p or- 
bitals plus one set of d orbitals on Si and Ge atoms, and 
two sets of s orbitals plus a set of p orbitals on H. To 
make the problem computationally treatable, most cal- 
culations were performed using single-C functions (SZ) 
for H as well as Si and Ge atoms at the three bottom 
slab layers, while DZP basis was kept on atoms at the 
remaining layers. Si and Ge atoms with SZ basis have 
one set of s and p orbitals, while H atoms have a single 
s orbital. Results obtained by using such a basis combi- 
nation are labeled with SZ-DZP. 

The unreconstructed bottom surfaces were saturated 
by hydrogen atoms making 1.50 A Si-H and 1.55 A Ge- 
H bonds. The opposite slab surfaces were set up ac- 
cording to specific surface structure models. These are 
dimer-adatom-stacking fault (DAS) based models, as well 
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as adatoms-based c(2 x 8), 2x2 and \/3 x \/3 recon- 
structions. The DAS models were initially proposed by 
Takayanagi et al. 13 for Si(lll)-7 x 7 and later extended 
for other DAS-family member structures (3 x 3, 5 x 5, 
7 x 7, 9 x 9, etc.). 14 Adatoms on c(2 x 8), 2 x 2, and 
^3 x V3 surfaces were placed at high-symmetry T4 ad- 
sorption sites. This is well established for adatoms on the 
Ge(lll)-c(2 x 8) reconstruction. 14 We also know that T 4 
adsorption sites are energetically more favorable than H 3 
and T\ in other cases. 15 The c(2 x 4) reconstruction, ob- 
served on quenched Si(lll) surfaces and on Ge/Si(lll) 
during growth was not considered. Our experimental re- 
sults indicate that it occupies a small fraction area be- 
tween different 2 x 2 or c(2 x 8) domain structures and 
therefore is considered a domain wall structure. 16 

A uniform real-space grid equivalent to a plane-wave 
cutoff of 200 Ry was used for Fourier transforming the 
density and potential terms. The geometry was opti- 
mized until all atomic forces became less than 0.01 eV/A, 
when the surface structure was considered to have at- 
tained equilibrium. All periodic slabs were separated 
by a 30 A thick vacuum layer. Under these conditions, 
converged calculations using bulk unit cells with a MP- 
S' 3 BZ-sampling scheme gave lattice constants of Si and 
Ge were asi = 5.420 A and clq c = 5.650 A. These are 
in reasonable agreement with the experimental values 
as; = 5.430 A and flQ C = 5.660 A, respectively. 



B. Convergence tests 

Here we describe the convergence tests to the condi- 
tions that have most impact on absolute surface forma- 
tion energies, starting with BZ sampling. To determine 
the optimum BZ sampling scheme for a particular slab 
lattice, we found the converged value of the Si (Ge) chem- 
ical potential from a bulk supercell employing a lattice 
identical to that of the slab (these geometries are here- 
after termed as bulk-slabs). To this end we progressively 
increased the density of the k-point grid until /1 changed 
by less than 5 /ieV/atom. Such a grid was considered to 
be appropriate for slab calculations with surface recon- 
structions (and vacuum gap). 

For Si and Ge bulk unit cells at zero strain, fj, con- 
verged at -107.844476 eV and -108.576450 eV, respec- 
tively, using a MP-8 3 grid of k-points. However, in order 
to reduce errors in relative surface energies, all calcula- 
tions used a common value of /1 obtained from a unitary 
bulk-slab, 3 layers thick and 1 x 1 in the (111) plane. 
With this geometry and MP-20 x 20 x 8, Si and Ge 
chemical potentials converged at fi = —107.844483 eV 
and n = —108.576456 eV, respectively, differing from the 
unit cell values by less than 7 /xeV/atom only. After care- 
fully looking at the convergence of fi using bulk-like slabs, 
we ended up with the following BZ sampling schemes for 
reconstructed surface slabs: MP-20 x 20 x 1 for 1 x 1, 
MP-6 x 6 x 1 for 3 x 3 DAS, MP-4 x 4 x 1 for 5 x 5 
DAS, MP-3 x 3 x 1 for 7 x 7 DAS, MP-2 x 2 x 1 for 
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FIG. 2. (Color online) (a) Si and Ge chemical potentials and 
(b) height factors, /h (see text) for Si(lll) and Ge(lll) in- 
terplanar spacings, calculated for various biaxial strain states 
eb in the (111) plane. 



9x9 DAS, MP-10 x 10 x 1 for 2x2, MP-8 x 2 x 1 for 
c(2 x 8) (rectangular surface cell) and MP-12 x 12 x 1 
for \/3 x \/3- The resulting k-point densities in recip- 
rocal space are approximately the same for all schemes. 
The estimated error in surface formation energies when 
compared to calculations with denser grids is well below 
0.1 meV/A 2 for all reconstructions. 

Figure 1(a) shows the calculated surface energy using 
DZP basis functions, for the ideal Si(lll)-1 x 1 (relaxed 
and unrelaxed) and reconstructed Si(lll)-3 x 3 DAS sur- 
faces as a function of the number of atomic layers in 
symmetric slabs. Beyond 12 layers the surface energy 
is converged within 0.01 me V/A 2 . These were therefore 
the number of layers used to calculate the reference sur- 
face energy 7i x i(e) from a Si(lll)-1 x 1 symmetric slab. 
Figure 1(b) depicts the Si(lll)-3 x 3 DAS surface en- 
ergy calculated using hydrogenated slabs with full DZP 
and mixed SZ-DZP basis as a function of biaxial surface 
strain. The replacement of DZP by mixed SZ-DZP basis 
causes a considerable decrease in the size of the Hamil- 
tonian, at the expense of a systematic downward shift of 
surface energies by about 0.5 meV/A 2 across all strain 
window of interest. Accordingly, we used DZP basis func- 
tions on Si (Ge) atoms belonging to the 3 upper layers, 
while cheaper SZ basis sets were assigned both to atoms 
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at the lowest 3 layers as well as to the hydrogen species. 
Although this step introduces the largest error in the cal- 
culated surface formation energies - about 0.5 me v/A 2 , 
it has little impact to the calculations of the surface strain 
response. 

Since most calculations were performed using hydro- 
genated slabs, we estimated the error produced by this 
approach. Usage of hydrogenated slabs requires less 
atomic layers, leading to a considerable reduction in com- 
putational effort in achieving self-consistent electron den- 
sities and energies. 8 We found that for hydrogenated 
slabs, 6 Si (Ge) layers were sufficient to provide highly 
accurate surface formation energies. The surface en- 
ergy difference between 6-layered hydrogenated and 12- 
layered symmetric Si(lll)-3 x 3 DAS surfaces is only 
0.03 meV/A 2 . Hence, the total number of atoms in hy- 
drogenated slabs were the following: 13 for lxl, 133 
for 3 x 3 DAS, 325 for 5 x 5 DAS, 641 for 7 x 7 DAS, 
1063 for 9 x 9 DAS, 53 for 2x2, 212 for c(2 x 8) and 40 
for v3 x \/3- For symmetric slabs the number of atoms 
was 30 for 1 x 1 and 248 for 3 x 3 DAS. The surface en- 
ergy for the 9x9 DAS reconstruction was calculated for 
zero strain only. A strain analysis for such a slab would 
require an incommensurate computational effort. 



C. Experimental procedure 

The experiments were performed in a ultra-high vac- 
uum system equipped with a STM (omicron). The 
12 x 3 x 0.4 mm 3 silicon samples were cut out from n-type 
silicon wafers with a resistivity of 0.5 Q, cm. The samples 
were heated by current injection and their temperature 
was controlled with help of an optical disappearing fila- 
ment pyrometer. 

The surface cleaning procedure was carried out at a 
pressure of 2 x 10~ 10 Torr. The Si(lll) samples were 
cleaned by degassing at least for 4 h at 600° C followed 
by flash annealing at 1250° C for a few seconds. We 
cooled the samples stepwise with 50° C per minute steps 
within the temperature range 900-350° C by reducing the 
heating current. Clean Si(7 7 10) surfaces were prepared 
following the procedure described by Kirakosian et al. 
Firstly, a sample was carefully degassed at 600° C for sev- 
eral hours. Then, after short flash annealing at 1250°C, 
the temperature of the sample was lowered to 1060°C 
during 30 s. After that the sample was quenched to 
830°C, kept at this temperature for 15 min, and finally 
cooled down to room temperature during 20 min. 

The source of germanium atoms consisted of pieces 
of Ge fastened to a W ribbon. The ribbon was heated 
upon current injection. Hereafter we specify the Ge de- 
posited coverage and thickness of Ge islands in bilayer 
units (1 BL Ge= 1.44 x 10 15 at/cm 2 ). These are equiv- 
alent to full Si (Ge) layers and must be distinguished 
from monolayers (ML), which are 1 ML = 1/2 BL for 
(111) surfaces of crystals with the diamond structure. 
This has been commonly adopted among experimental- 
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FIG. 3. (Color online) Surface formation energies calculated 
for various biaxial strain states eb and surface reconstructions 
for (a) Si(lll) and (b) Ge(lll) surfaces. 

ists working with Ge/Si(lll). The Ge deposition rate 
in our experiments was 0.001-0.01 BL/min. The Ge flux 
was calibrated by measuring the total volume of Ge is- 
lands formed on Si(lll) far from step edges and domain 
boundaries of surface reconstructions. 18 

STM images were recorded at room temperature in the 
constant-current mode using an electrochemically etched 
tungsten tip. The freely available wsxm software was 
used to analyze and process the STM images. 19 



III. RESULTS AND DISCUSSION 

A. Surface energies and stresses of Si(lll) and 
Ge(lll) 

With the exception of the c(2 x 8) surface structure, 
for which a small anisotropy on its surface stress tensor is 
expected, all other surface reconstructions in this study 
have C3 symmetry. For such surfaces the surface stress 
tensor is isotropic: a xx = a yy = a ICC and a xy = a yx = 0, 
where the [111] direction is assumed to be parallel to 
the z-axis. Hence, we applied an isotropic biaxial strain 

€xx — €yy 

= £b, and from the surface formation energy 
7rec( e t>) we found the surface stress cr r cc(eb = 0). For 
surfaces having C3 symmetry the dependence of 7 roc on 
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uniaxial strain e u can be achieved from 7 re c(eb) by using 
the relation, 

1 + e u = (1 + e b ) 2 . (4) 

Figure 2(a) shows the calculated chemical potential 
fi as a function of (111) biaxial strain, £b, for bulk Si 
and Ge. From elasticity it follows that a (lll)-biaxially 
strained cubic solid suffers an opposite strain along [111]. 
Such effect has to be accounted for in strained surface cal- 
culations, and so one should let the surfaces relax towards 
the vacuum. This effect was considered by using appro- 
priately strained lxl bulk-slabs with several heights, 
h, related to the strain- free height />o by a height factor, 
fh = h/h,Q. We determined their equilibrium heights [h 
values that minimized the energy) , which were then used 
to obtain the chemical potentials of Si (Ge) for a specific 
strain state. Fig. 2(b) depicts calculated height factors 
for Si and Ge 1 x 1 bulk-slabs. It shows how the equi- 
librium distance between (111) layers in bulk depends on 
the applied biaxial strain. A direct comparison between 
Si- and Ge-related plots clearly shows that in agreement 
with the experimental data, Si has greater mechanical 
stiffness than Ge. Both Fig. 2(a) and Fig. 2(b) show that 
the plots are slightly asymmetric and hence the strain in- 
terval spans the anharmonic regime. 

Figures 3(a) and 3(b) show the calculated surface en- 
ergies 7roc(eb) of Si(lll) and Ge(lll), respectively, for 
various experimentally observed reconstructions. These 
were calculated using combined SZ-DZP basis for e b val- 
ues between —0.04 to 0.04, and shifted by their respec- 
tive systematic errors with respect to the more complete 
DZP basis calculations at e b = 0. Calculations of the un- 
reconstructed Si(lll)-1 x 1 are not shown, for this surface 
has never been observed experimentally. Note that the 
Si(lll)-1 x 1 structure observed above the 7 x 7 — > 1 x 1 
phase transition temperature is not related to the static 
lxl structure reported here. It is actually caused by 
a thermally-average motion of atoms. The surface en- 
ergy of the ideal (relaxed) lxl surface monotonically 
decreases from 127.1 to 102.0 meV/A 2 for Si(lll), and 
from 93.5 to 77.6 meV/A 2 for Ge(lll) when increasing 
the surface strain from —0.04 up to 0.04 (compressive sur- 
face stress). Surface energies of all other reconstructed 
surfaces show an opposite trend. They increase with 
strain (tensile surface stress) as seen in Fig. 3. According 
to elasticity theory of solids, a non-vanishing derivative 
of the energy with respect to strain at e = indicates 
an intrinsic driving force for explosion or implosion of 
the crystal. This is not the case for surfaces since they 
couple to the underlying bulk, and therefore are not free 
to expand or contract infinitely in response to a surface 
stress. 

Looking at the surface energy plots for Si(lll) in 
Fig. 3(a), they can clearly be grouped into two sets. The 
first set includes 3 x 3, 5 x 5 and 7x7 DAS reconstruc- 
tions, while the second set includes c(2 x 8), 2 x 2, and 
Vox v3 structures. As explained below, such distinction 



TABLE I. Absolute surface energies and stresses (meV/A 2 ) 
for various Si(lll) and Ge(lll) surface reconstructions cal- 
culated using DZP basis at zero strain. Previously reported 
data are also shown for comparison (footnoted values). 

Energy, 7 rcc Stress, cr rec 
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Ge 
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116.95 

1 1 o C'ra 
llo.O / 


86.45 

OZ.O ( 








112.33 


85.53 


-42.94 


-27.09 


lxl, relaxed 


lUs.OS 

113. 99 b 


ol.lz 


QO Q1 b 




6X6 JJAo 


f\A Q £ 

94.80 
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210. U I 
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5x5 DAS 


92.27 


68.83 


201.72 


141.28 


7x7 DAS 


92.42 


68.73 


191.80 


135.23 


84.86 a 


63.65 a 






c(2 x 8) 


95.00 
87.98 a 


68.62 
63.02 a 






2x2 


95.32 
105.34 b 


68.79 


144.38 
141. 51 b 


95.65 


x 


95.70 


71.09 


141.23 


98.40 



a Prom Ref. 8. 

b From Ref. 6, assuming agi = 5.420 A. 



stems from different building blocks that make their re- 
spective reconstructions. Analogous sets are observed in 
Fig. 3(b) for Ge(lll) as well. The only exception being 
that of the \/3 x \/3 surface energy data (adatom-based 
structure with the highest density of adatoms), which 
is somewhat less stable than other members of the sec- 
ond group. Comparing 7 rcc (£b) curves from Figs. 3(a) 
and 3(b), one readily conclude that the formation of sim- 
ple adatom-based Ge(lll) reconstructions compete with 
DAS-based counterparts near the region of zero strain, 
and that each structure type is prominent depending on 
the strain sign. On the other hand, Si(lll) DAS-based 
structures are the most stable across most of the strain 
interval. This indicates that additional Si reconstruction 
elements like dimers and stacking faults are particularly 
important in the stabilization of the Si(lll) surface. 

In Table I we report absolute surface energies (DZP 
basis) and stresses (DZP-corrected SZ-DZP calculations) 
for various Si(lll) and Ge(lll) surface reconstructions 
at zero strain. Surface energy curves, 7 rec (e)j were cor- 
rected to match DZP results at e — 0. Surface stress 
values were taken from the first-order variation of 7 rec (e) 
at equilibrium, 3 



<7 r ec=(jrec + ^J 6=o , (5) 

where e u is the uniaxial strain (that relates to e b through 



6 



TABLE II. Absolute surface energies 7 rcc at zero-strain 
(me V/A 2 ) for DAS-based reconstructions of Si(lll) and 
Ge(lll) surfaces calculated using the combined SZ-DZP ba- 
sis. 



Reconstruction 


Si 


Ge 


3x3 DAS 


94.37 


70.11 


5x5 DAS 


91.88 


68.57 


7x7 DAS 


92.06 


68.49 


9x9 DAS 


92.48 


68.58 



Eq. 4), and 7 rcc ( e u) wa s approximated to a quadratic 
function that was fitted to the first-principles data. 

Table I shows that unstrained Ge(lll) with c(2 x 8) 
reconstruction is correctly predicted to be the ground 
state. Indeed, the experiments show that clean (111) 
surfaces of Ge samples have a c(2 x 8) structure. 14 

According to our calculations, strain-free Si(lll) sur- 
faces with 5x5 and 7x7 DAS reconstructions are the 
most stable and they are almost degenerate. There is 
however a small preference for the 5x5 reconstruction 
by about 0.15 meV/A 2 , which is just beyond the esti- 
mated error of 0.1 meV/A 2 . The obtained ordering is 
at variance with the experiments. Clean Si(lll) samples 
obtained after a high temperature annealing and a slow 
cooling stage down to room temperature, normally show 
a 7 x 7 DAS reconstruction. It should be noted though, 
that there is indirect experimental evidence that the sur- 
face energy of Si(lll)-5 x 5 is rather low and close to 
that of the 7x7 surface. Indeed, the 5x5 reconstruc- 
tion has often been observed close to surface defects such 
as steps. 20 Further evidence was found during the phase 
transition from the metastable 2x1 reconstruction to 
7x7 DAS, which proceeds through the formation of an 
intermediate 5x5 DAS surface structure. 21 

The energies of 3 x 3, 5x5 and 7x7 DAS recon- 
structed Si(lll) surfaces were also calculated using the 
DFT method in Ref. 22. The authors found that the 7x 7 
DAS reconstruction was the ground state, only 2 meV/A 2 
and 3 me v/A 2 more stable than 5x5 DAS and 3x3 DAS 
surface reconstructions, respectively. However, these cal- 
culations were shown to suffer from poor BZ sampling, 
that once adequately corrected, resulted in a 5 x 5 DAS 
ground state 5 me v/A 2 more stable than the 7x7 DAS 
reconstruction. 23 Still, it is not clear why despite the 7x7 
being the most commonly observed surface reconstruc- 
tion in Si(lll) after thermal treatment, the calculated 
Si(lll)-5 x 5 surface energy is actually predicted to be 
the ground state. Among the possible reasons for such 
disagreement one could have an insufficient number of 
layers considered in the slab, a poor treatment of the 
exchange and correlation energy, the neglect of entropy 
effects or the need for a better description of the electron 
density and potential. 

To have an idea of the effect of the basis-choice to the 
surface energy, Table II lists absolute surface energies of 
various DAS reconstructions on Si(lll) and Ge(lll) cal- 



TABLE III. Calculated Si(lll) surface stresses <r (meV/A 2 ) 
for various experimentally observed reconstructions, com- 
pared to those of conceived model structures. 



Observed 




Test models 




Reconstruction 


o 


Reconstruction 


a 


1 x 1 


— 42.94 


1 x 1 SF 


— 1.06 


2x2 


144.38 


2 x 2 SF 


163.85 


X 


141.23 


V3x^3SF 


163.81 


3x3 DAS 


215.07 


3 x 3 DS 


174.43 


5x5 DAS 


201.72 


5 x 5 DS 


138.98 


7x7 DAS 


191.80 


7 x 7 DS 


113.73 



culated with SZ-DZP basis at zero strain. From a com- 
parison with Tab. I we conclude that the energy difference 
between Si(lll)-7 x 7 and Si(lll)-5 x 5 reconstructions 
decreases by less than 0.5 meV/A 2 upon increasing the 
basis size. This result may direct us on future studies. 

The magnitude and sign of the surface stress is di- 
rectly related to the details of chemical bonding at the 
surface. It is possible to obtain a qualitative estimate 
of the impact of specific reconstruction elements to the 
overall surface stress. This is done by conceiving and in- 
specting test models that include elemental structures of 
interest, even if the resulting surfaces are not observed 
experimentally. To this end the following test models 
were produced (i) 1 x 1 SF, 2 x 2 SF and ^3 x ^3 SF re- 
constructions, where stacking faults (SF) are introduced 
between the first and second surface layers, and (ii) 3x3, 
5x5 and 7x7 dimers stacking fault (DS) structures, that 
basically are DAS counterparts after removing adatoms. 

The comparison between surface stress calculations for 
observed and test model reconstructions is summarized 
in Tab. III. Several conclusions can be drawn from the 
data - firstly, comparing the stresses of 1 x 1, 2x2 
and y/3 x v3 Si(lll) surfaces with their counterpart val- 
ues from SF structures, one may conclude that stacking 
faults produce a tensile contribution to the surface stress. 
Secondly, given that both DAS and DS structures have 
dimers located at the boundary of half unit cells, when 
their size is increased from 3 x 3 to 5 x 5 and then to 7 x 7 
the relative contribution of dimers to the surface stress 
obviously decreases. This is accompanied by a decrease in 
the surface stress, meaning that the dimers in DAS recon- 
structions also produce a tensile contribution. Thirdly, 
when adatoms are added to the Si(lll) truncated lxl 
surface to produce 2x2 and y/3 x reconstructions, 
they produce a strong tensile stress. Conversely, when 
adatoms are removed from surfaces having DAS recon- 
structions, the tensile strain decreases. 

In summary, all reconstruction elements, i.e. dimers, 
adatoms and stacking faults produce a tensile contribu- 
tion to the surface stress. A tensile contribution to the 
surface stress is also expected from the reconstruction ele- 
ments on (111) surface of germanium. The origin of com- 
pressive stress in the truncated Si(lll) and Ge(lll) with 
lxl periodicity lies at the Coulomb repulsion among the 
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FIG. 4. (Color online) STM images obtained during different formation stages of a germanium wetting layer on Si(lll) at 
T = 400° C (quasi-equilibrium growth conditions), (a) Ge island formed during the initial stage exhibiting a 7 x 7 reconstruction 
with y/~S/H — 11.31. At this stage, the Ge total coverage was 0.6 BL. (b) On rare occasions one can observe Ge islands with a 
y/3 x y/i reconstructed top. Here y/S/H = 24.33 and the Ge total coverage was 0.4 BL. (c) During lateral growth, top surface 
reconstructions transform to 5 x 5. The dashed line highlights the domain boundary between 7x7 (left) and 5x5 (right) 
surfaces. The Ge total coverage was 0.8 BL and for this island y/S/H = 31.68. All islands in the figures are 3-BL thick. 



two-dimensional array of surface dangling bonds. 3 We 
also note from Fig. 3 that the expected structure change 
is c(2 x 8) — > 7x7— > 5x5 with increasing compres- 
sive strain for both Si(lll) and Ge(lll) surfaces. This 
is because compression relieves the intrinsic tensile stress 
so that reconstructions with bigger tensile stress are the 
most stabilized. 



B. Experimental results and comparison with 
theory 

We begin by looking at surface structure transfor- 
mations that take place during the early stages of 
Ge/Si(lll) growth. It is useful to define the parameter 
y/S/H to characterize the shape of growing Ge islands, 
where S is the area of their top terraces (y/S is a charac- 
teristic lateral size of the terrace) , and H is the height of 
the island. Thus, y/S/H relates to the ability of the is- 
land to relax elastically, that is, higher values correspond 
to islands with a highly strained Ge/Si(lll) upper layers, 
while lower values correspond to a more relaxed one. 

The growth of Ge on Si(lll) proceeds according to 
the SK growth mode, consisting in the formation of a 
few monolayers thick wetting layer during a first stage, 
followed by the growth of individual 3D islands as fur- 
ther deposition proceeds. 3 ' 4 The thickness of the wet- 
ting layer depends on the rate of Ge deposition and 
temperature. 24 When growing at quasi-equilibrium con- 
ditions, as it was the case for this work (slow deposi- 
tion rate 10 -3 -10 -2 BL/min and high temperature, i.e. 
300-500° C), the formation of the Ge wetting layer takes 
place upon growth and coalescence of neighboring 3-BL 
thick Ge islands. If a higher deposition rate and/or a 
lower temperature are used, the formation of Ge islands 
with such thickness becomes kinetically limited, and two- 
bilayer or even single-bilayer thick islands may form. 



Figures 4(a)-(c) show STM images of typical 3-BL 
thick Ge islands grown at T = 400° C at various stages of 
the wetting layer formation. The STM image of a typ- 
ical island formed at the very beginning of Ge/Si(lll) 
growth is shown in Fig. 4(a). The strained Ge islands 
formed at this stage exhibited the 7x7 reconstruction 
pattern of the underlying Si substrate surface. Rarely 
at this stage one could find Ge islands with a y/3 x y/3 
reconstructed top terrace as it is depicted in Fig. 4(b). 
Following further growth, the islands increase their lat- 
eral size (although maintaining their heights), and upon 
reaching some critical size, the 5x5 reconstruction pat- 
tern appears on top terraces [see Fig. 4(c)]. Finally, the 
Ge islands merge together to form a continuous film (the 
so called the wetting layer) , fully 5x5 reconstructed with 
no traces of 7 x 7. 

Figures 5(a)-(c) show STM images of typical 3D Ge is- 
lands grown at T — 400° C on top of the Ge wetting layer 
during various stages. Small 3D Ge islands formed early 
have well developed side facets and tiny top terraces. 25 
One of such 3D islands is shown in Fig. 5(a), with its 
top terrace exhibiting c(2 x 8) and 2x2 domain recon- 
structions. Upon further growth, these 3D Ge islands 
increase in size, both laterally and vertically. Some of 
them exhibit a 7 x 7 reconstruction on their top surfaces 
[see Fig. 5(b)], while others show a c(2 x 8) and 2x2 
structures as shown in Fig. 5(c). Note that the obser- 
vation of c(2 x 8) reconstructed top layers indicates that 
they reached a relaxed state. This is the strain-free stable 
structure for clean (111) surfaces in pure Ge. 14 

The calculated surface phase diagram shown in 
Fig. 3(b), accounts rather well for the observed sur- 
face structure transformation sequence that was just de- 
scribed. The Ge islands that are formed before com- 
pletion of the wetting layer are compressively strained 
due to lattice mismatch (although some strain must be 
relieved at their edges). According to Fig. 3(b) the ef- 




FIG. 5. (Color online) STM images obtained during different stages of Ge growth on a Ge/Si(lll) wetting layer at T = 400°C 
at quasi equilibrium conditions. This layer was invariably 5x5 reconstructed, (a) Ge islands formed during the early stages 
show small top terraces with c(2 x 8) and 2x2 reconstructions. The Ge total coverage here was 2.5 BL and for this particular 
island ^fS/H = 0.71. (b), (c) Upon further growth, and depending on the ^fS/H ratio, one observe both 7x7 and c(2 x 8) 
reconstructed terraces. The islands depicted have V~S/H = 13.15 and \f~S/H = 5.15, respectively, and the Ge total coverage 
was 3.0 BL in both cases. 



fective built-in strain in such islands should be within 
—0.01 < eb < 0.0. In this strain range the 7 x 7 re- 
construction is the most favorable state. As the islands 
increase in size, their ability to relax laterally decreases. 
Thus, wider Ge islands and particularly the continuous 
wetting layer, accumulate more compressive strain and 
that induces the 7 x 7 — > 5 x 5 phase transition in agree- 
ment with Fig. 3(b) for -0.04 < e b < -0.01. 

When it comes to the formation of 3D Ge islands on 
top of the wetting layer, it is important to follow their 
V S / H ratio. Initially, only islands with tiny top terraces 
and very small VS/H values can be found. Due to their 
small lateral dimensions, these islands effectively release 
the compressive strain at their top layers. According to 
Fig. 3(b), their c(2 x 8)-reconstructed terraces indicate 
that they are mostly relaxed or eventually could accumu- 
late some tensile strain. Larger 3D islands show cither 
c(2 x 8) or 7 x 7 reconstruction, and again, we find a clear 
correspondence between \/~S/H ratios and the observed 
terrace reconstructions: islands with lower ^/S/H show 
relaxed c(2 x 8) terraces while islands with higher / H 
exhibit compressive 7x7 terraces. 

Our results indicate that the 2x2 reconstruction is 
metastable across the studied strain domain. Its forma- 
tion during Ge/Si(lll) growth, must be then driven by 
kinetics (not considered in this work) . This structure was 
however frequently observed along with the c(2 x 8) re- 
construction, and this is in line with the small calculated 
formation energy difference between these surfaces. The 
Ge(lll)--\/3x y/3 reconstruction with adatoms at T4 posi- 
tions is strongly unfavorable across the whole —0.04-0.04 
strain range. The surface energy was therefore recalcu- 
lated assuming H3 high-symmetric positions for adatoms, 
although resulting in an even less stable structure. It is 
possible that the y3 x -\/3 structure is not made of Ge 
species only, and some Si intermixing could play an im- 
portant role on its stabilization. 



Another interesting system to observe the influence of 
strain on surface reconstructions is a stepped surface. 
The high-index regular stepped Si(7 7 10) surface, ideally 
consists of (111) terraces comprising half of 7x 7 unit cells 
in width and steps with height of three (111) interpla- 
nar spacings. 26 The regularity of this surface makes it an 
ideal substrate for growing nanodots and nanowires. 27 It 
should be noted that its orientation and atomic structure 
was the source of much controversy as it was formerly 
identified with a Si(557) surface. 17,28 Presently, the most 
consensual atomic model and the one that agrees better 
with the available experimental data is that reported in 
Rcf. 26. 

Whereas Si(lll) flat samples demonstrate surfaces al- 
most fully 7 x 7-reconstructed, the (111) terraces of the 
Si(7 7 10) stepped surface sometimes show 5x5 struc- 
ture patches. 20 This is depicted in Fig. 6(a). According 
to Fig. 3(a), the observation of a 5 x 5 structure means 
that the surface layers at the (111) terraces suffer a com- 
pressive strain. Indeed, unlike the flat (111) samples, the 
(111) terraces on stepped surfaces have a nearby free bor- 
der at the step edges. The Si(lll)-7 x 7 surface exhibits 
tensile stress (Tab. I), therefore in the presence of a free 
border, the surface area tends to shrink. This effect leads 
to a shift of the strain state as depicted in Fig. 3(a) to- 
wards the compressive side (leftwards) , making the 5x5 
surface structure the most stable. This reasoning is, how- 
ever, qualitative as it overlooks the structure of the step 
edges. These should differ for 7 x 7 and 5x5 recon- 
structed (111) terraces. The 5x5 reconstruction of a Ge 
layer deposited on Si(7 7 10) can be explained based on a 
similar reasoning [see Fig. 6(b)]. 29 

Finally we note that theory predicts that Si(lll) sur- 
faces with considerable tensile strain (above 0.03) should 
show a c(2 x 8) reconstruction [cf. Fig. 3(a)]. Such 
strain can be produced when growing Si on Ge(lll) sur- 
faces and this should be verified in future experimental 
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FIG. 6. (Color online) STM images of a stepped Si(7 710) 
surface with (111) terraces, (a) Clean Si(7 710) surface with 
triple steps. All terraces show a 7 x 7 reconstruction except 
the one on the left, which has 5x5. (b) Si(7 710) surface 
after adsorption of 0.3 BL Ge at T = 600°C. Triple steps are 
split into single and double steps. All terraces exhibit a 5 x 5 
reconstruction. 



wide set of experimentally observed surface reconstruc- 
tions. The calculations were carried out within the lo- 
cal density functional framework, by applying several bi- 
axial strains to vacuum-separated slabs with appropriate 
surface reconstructions. Using these data, the surface 
stresses of Si(lll) and Ge(lll) having different recon- 
structions, were calculated at the strain-free state. The 
results indicate that dimers, adatoms and stacking faults 
all produce a tensile contribution to the surface stress. 

The formation energy of the Si(lll) unstrained surface 
with a 9 x 9 DAS reconstruction was also calculated. It 
is shown that the energy of this surface is indeed higher 
than corresponding energies of surfaces having 5x5 and 
7x7 DAS structures. However, the Si(lll)-7 x 7 surface 
is found to be slightly higher in energy than the Si(lll)- 
5 x 5, in apparent contradiction with the observations. 
This inconsistency calls for further work on this problem. 

Strain phase diagrams for Si(lll) and Ge(lll) sur- 
faces were constructed, showing a direct relation between 
the equilibrium surface structure and the applied strain. 
Comparison of the theoretical results with STM data 
shows that the constructed phase diagrams can explain 
the surface structure transformations observed during 
growth of strained Ge/Si(lll), as well as the resulting 
structures found on the stepped Si(7 710). Therefore, 
the observation of surface reconstructions may be used 
as a probe to the local strain, provided that we have a 
surface-strain phase diagram like the one presented in 
Figure 3. Finally, theory predicts that Si(lll) should 
adopt a c(2 x 8) reconstruction in the presence of a ten- 
sile strain approximately above 0.03. 
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